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Abstract 

Multicolor photometric observations of a planetary transit in the system HD 209458 
are analyzed. The observations, made in the Stromgren photometric system, allowed 
a recalculation of the basic physical properties of the star-planet system. This in- 
cludes derivation of linear limb-darkening values of HD 209458, which is the first 
time that a limb-darkening sequence has observationally been determined for a 
star other than the Sun. As the derived physical properties depend on assumptions 
that are currently known with limited precision only, scaling relations between de- 
rived parameters and assumptions are given. The observed limb-darkening is in 
good agreement with theoretical predictions from evolutionary stellar models com- 
bined with ATLAS model atmospheres, verifying these models for the temperature 
(T eff 6000ET), surface gravity (log g « 4.3) and mass (« 1.2M ) of HD 209458. 
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1 Introduction 



The discovery of the first extrasolar planet transiting its central star, HD 209458 
(Charbonneau et al., 2000; Henry et al., 2000) opened a new way to study 
these objects through the analysis of the light curve that is generated during 
a transit. Similarly, details about the central star may be revealed. Here we 
present observations of a transit of HD 209458 that were taken simultane- 
ously in the four Stromgren colors. The primary motivation was a study of 
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the specific intensity distribution of the stellar surface of HD 209458 for com- 
parison with existing models of stellar atmospheres through the derivation of 
limb-darkening coefficients. A transit light curve provides in its central part - 
between second and third contact - a direct measure of the specific intensity 
behind the planet. The specific intensity behind the ingress and egress zones 
is also accessible, through models of transit light curves. A planetary transit 
profile gives significant advantages against efforts to derive intensities across 
the stellar disc from light curves of eclipsing binaries: in transits, a small dark 
disk passes in front of the star, whereas eclipsing binaries cause ambiguities 
from the presence of two luminous sources, frequently of similar size range, 
and often aggravated by light reflection between the two components, or by 
gravitational darkening in the case of fast rotators. 

The major difference between the data presented here and the BVRIZ multi- 
color transit observations by Jha et al. (2000) is the different wavelength span, 
with our data taken in the Stromgren photometric system. Also, our data were 
taken in all colors simultaneously, whereas Jha et al. used a CCD camera cy- 
cling sequentially through the V, R and / filters, with data in further colors 
obtained at telescopes in different locations. In the following analysis we then 
show, that transit light curves do not only allow the derivation of parameters 
of the extrasolar planet, but impose severe constraints on the mass-to-radius 
relationship of the central star as well. Furthermore, derived parameters of the 
star-planet system are presented in a way that is independent of assumptions 
based on approximately known values, such as the stellar mass. Finally, linear 
limb-darkening coefficients are given as a function of the stellar radius and 
mass of HD 209458 and compared to stellar and atmosphere models. This is 
the first time that a color sequence of limb-darkening coefficients for a single 
star besides the Sun has been derived from observations. 



2 Data acquisition 

The observations of HD 209458 reported here were performed in the night of 
25-26 July 2000 with the 0.9m telescope of the Observatorio Sierra Nevada, lo- 
cated near Granada, Spain at an altitude of 2900m. The telescope is equipped 
with a 4 channel photometer, allowing simultaneous observations in all colors 
of the u, v, b, y Stromgren photometric system. This instrument is similar to 
that of the Stromgren Automatic Telescope at La Silla (Gronbech et al., 1976). 
The observing cycle consisted of a 30 sec sky-exposure, 20 sec exposures of 
standard stars SAO 107675 and SAO 107561, and a 40 sec exposure of the 
object, HD 209458. During the ingress and egress of the transit, for about 40 
minutes each time, only HD 209458 was observed with repeated 40 second 
exposures. During the flat bottom part of the transit, the second comparison 
star was left out. These modifications served to increase the time resolution, 
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Fig. 1. u, v, b, and y light curves that were used as the basis for the analysis (crosses). 
They are smoothed by a triangular kernel with a width of 5 data points. The solid 
lines are model light curves for a star-planet system with the parameters given in 
Table 3. 
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and were justifiable since the night was photometric, with last quarter moon, 
and the mid-transit occurred within a few minutes of the meridian crossing 
of HD 209458, leading to equal airmasses at ingress and egress. Data reduc- 
tion was performed at the telescope Macintosh control computer, with the 
derivation of an extinction function based on the comparison stars. Differ- 
ential photometry was then performed with respect to the first comparison 
star. Resulting lightcurves are available in electronic Table 1. Magnitude dif- 
ferences between both comparison stars indicate a photometric precision of 
about 4 mmag for the u and y bands and 2 mmag for the v and b ones per 
single measurement, with the better results in the two latter bands simply 
being due to higher photon count rates. After removal of poorer data at the 
begin of the observations, the transit light curves used in the analysis contain 
129 data points in all colors, with 29 off-transit points, which were used to set 
the magnitude zero-point 2 . In the off-transit part, the rms noise is: 4.2mmag 
for u, 2.2mmag for v, 1.9mmag for b and 3.4mmag for y. For the subsequent 
analysis, the light curves were smoothed with a triangular kernel of 5 elements 
width, obtained through two successive boxcar smoothings with a 3 element 
rectangular kernel. The off-transit rms of the resulting light curves, shown in 
Fig. 1, is then: u: 1.8mmag, v: 0.72mmag, b: 0.89mmag, y: 2.5mmag. Some 
parts of our analysis are based on a combined v + b light curve, for which a 
lower rms of 0.65mmag was attained. Similar observations as reported here 
were performed in the night of July 18-19, 2000 as well, but increasingly hazy 
weather conditions and a nearly full moon significantly degraded these obser- 
vations, leading us to concentrate on the data from 26 July. 



3 Analysis of the transit 

3. 1 Setup of model fits to the transit 

A convenient set of 'primary transit parameters', from which other frequently 
used parameters (like the planet's orbital inclination % or the orbital half axis) 
can easily be derived, may be selected as follows: 

P: orbital period 

T c : epoch of the planetary transit 

M*: stellar mass 

i?* : stellar radius 

Rpi'. effective planetary radius, which is the radius of an opaque disk occulting 

2 For the y filter, shorter term variations not well described by the extinction func- 
tion required offsets from the off-transit zero-point found through adjustment of the 
best fitting planet size (see Fig. 6) to that of the other colors 
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the same amount of light as the planet. 

b: latitude of the transit across the stellar disk in stellar surface coordinates 
w: linear limb darkening coefficient. 

w is representative for some descriptor of the intensity distribution across the 
stellar disk. Here, we will use the common linear limb-darkening expression, 
defined through I (cos 6) = 1(0) [1 — w(l — cos0)], where 9 is the angle between 
the stellar surface and the line of sight. The use of a linear limb-darkening 
expression for HD 209458 is justified, since it gives a good description of surface 
intensities for stars with effective temperatures between 5000 and 6000 K, and 
T eff = 6000 K for HD 209458 (Mazeh et al, 2000). 

Among the primary parameters, only the orbital period and the epoch are 
known with high precision (Soderhjelm, 1999; Castellano et al., 2000; Robichon 
and Arenou, 2000). For this work, a period of 3.52474 days and an epoch 
of HJD 2451659.93675 from HST observations by Brown et al. (2001) was 
used, and - upon excellent agreement with our data - no attempt for a new 
derivation has been undertaken. Furthermore, estimates with a precision of 
about 10% exist for and M*, with = 1.2R Q and = l.li2 (Mazeh 
et al., 2000, and references therein). 

A boundary condition that allows a reduction in the number of free parameters 
is given by the total length of the transit, which is the time T i4 between 
the first and the last - fourth - contact of star and transiting planet. The 
observations by Brown et al. (2001) give a value for T 14 of 0.12795 days, 
which is in excellent agreement with our data. T 14 is given by the following 
expression, whose derivation is shown in Fig. 2: 



where the planet's orbital velocity t> or b depends in circular orbits 3 on M* and 
P as f or b oc (M^/P) 1 ^. With T i4 known, Eq.(l) can be solved for cos b as a 
function of i?*, M* and R p \: 



which corresponds to physically meaningful solutions only for < cos b < 1. 



3 For the case of HD 209458b - a planet that is on a very close orbit around its 
central star - the actual T 14 will be, in the worst case, 1.0024 times longer than 
derived by Eq.(l) 




(1) 
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Fig. 2. The total transit duration T14 = <Str/v rb depends on the lateral length of the 
transit S tr as follows: (St r /2) 2 + h 2 = (R p \ + R*) 2 , with R p \ and i?* being the radii 
of the planet and of the star. If h < i?*, then h may be expressed in terms of the 
latitude b of the transit across the surface of the star: h = R* sin b, and consequently, 
h 2 = R 2 (l - cos 2 b). Solving this for S tl leads then to Eq.(l). 

For a given combination of the parameters cos 6, i?*, i? p i, an estimate for the 
limb-darkening coefficient w can then be obtained as follows: 
The relative brightness change AL/L at the center of a transit, assuming the 
linear limb-darkening expression, is given by: 



which can be derived from the linear limb-darkening equation, with 6 = b at 
the center of a transit. With an estimate of AL/L taken from the observed 
light curve, the above equation can be solved for w: 

w = 3 j(R pl /R*) 2 - AL/L] 

3(R pl /R*) 2 (l - cos b) - AL/L { ' 



Fits to the light curve were performed using the transit modeling code UTM, 
which calculates light curves following a set of input parameters, and the 
fitting routine 'UFIT' 4 , which wraps around UTM and provides UTM with 

UTM (Universal Transit Modeler) and UFIT (Universal Fit- 
ter) are written in IDL v. 5 code and are available from HJD at 
http://www.iac.es/galeria/hdeeg/idljians_lib/utm/. UTM can be used 
in stand-alone mode to produce transit light curves, or within the fitting routine 
UFIT. UTM calculates the light curves of transits from pixelized representations 
of projections of the stellar surface and the transiting object. For both UFIT and 
UTM, setup files are used to completely describe a transit configuration, and allow 
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Fig. 3. The contours give limb darkening values 'w' for the combined v + b color, as 
a function of the boundary parameters R p \/R* and i?*, for a fixed stellar mass of 
1.1 Mq (or Rpi/R* and R'^ for any stellar mass, see Sect. 3.2). The limited vertical 
lengths of the contours is given by the range in which solutions where calculated. 
The dashed line gives the trace of the best fit for each value of R'^ (see also Fig. 4). 
Small changes in the planet size (up-down) around the best fit have little effect on 
w. 

input parameters following a Levenberg-Marquardt fitting algorithm. UFIT 
then stepped through a 3 dimensional grid of boundary values for M*, R* 
and Rp\/ R*, performing fits of the only free parameter w at each step. The 
boundary values were varied within the following limits: 

0.95M Q < < 1.35M , with increments of O.O25M 
O.95i? < R* < 1.35^0, with increments of O.O25i? 
0.11 < Rpi/R* < 0.145, with increments of 0.00125 

At each step, a new value for cos b was calculated from Eq.(2). If a physically 
meaningful solution for cos b was found, an initial value for w was calculated 
from Eq.(4). Only if this value was in the range of 0.2 < w^a < 1, an actual fit 
to the light curve was performed. Since w does not appear in Eq.(l), results 
from the fits will always comply with the boundary condition of the total 
transit duration Ti 4 . 

the representation of systems with several stars and planets, moons or rings. UFIT 
can perform fits to any of the input parameters of UTM. 
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Fig. 4. Confidence regions (1, 3 and 5 cr) of the best fit to w, within the same 
parameter space as the previous figure. The dashed line traces the best fitting value 
of R p i/R* for a given R'^. Dotted lines indicate the latitude (in degrees) of the 
transit, from Eq.(2). 

The result from these 'scanning fits' are two three-dimensional arrays with 
the axes along M*, R* and Rp\/R*, one containing values for the fitted limb- 
darkening coefficient w and the other one the variances between fit and data. 
A 'cut' through such an array of w values, from fits to the combined v + b 
light curve, is shown in Fig. 3, and the corresponding variances, converted 
to confidence regions, are shown in Fig. 4. All confidences in this paper are 
calculated for 129 data points and 4 free parameters. 

3.2 Dependence on stellar mass and radius 

In the following, we will use the light curve with the lowest noise -the combined 
one from v and b niters- for the derivation of all parameters, except for the 
limb-darkening coefficients, which are wavelength dependent. In Fig. 5, the 
3-dimensional array that resulted from the scanning fit is flattened into 2 
dimensions (M*, R*) by selecting those values of R p \/ R* which gave the best 
fit at each point M*,R*. The best solution appears in a narrow band, which 
can be approximated by 



= 0.34M* + 0.825(±0.06), 



(5) 
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Fig. 5. 1, 3, and 5 a confidence regions of the best fit on the combined v and b light 
curve, as a function of the mass and radius of the central star. Best fits lie along a 
band in (M*, R*), also given by Eq.(5). The 5 points with error bars are those radii 
where limb-darkening coefficients from transit data and from ATLAS atmosphere 
models coincide (see Section 4 and Fig. 8). 

with the errors indicating a 3a confidence limit. This relationship follows lines 

-| /o 

of constant transit latitude, which is dominated by cos b oc v or h/R* oc Mj / R* 
for R p i « R* (see Eq.2). It should be noted that this relationship has no 
physical meaning in the sense of compatibility with stellar models. It is simply 
a result from the transit light curve and the T 14 boundary condition, giving 
the range of possible radii as a function of the stellar mass that are compatible 
with transit observations. Eq.(5) allows us to express most results given in the 
further sections as functions of a 'normalized radius' R'^, which is related to 
the actual stellar mass and radius (currently known with large errors only) as 
follows: 

K = R* + 0.34(1.1 -M*) (6) 

For M* = 1.1, normalized and actual radii are equal. It should be noted that 
the stellar parameters given by Mazeh et al. (M* = l.lM ,i2* = 1.2i2 ) lie 
well within the la confidence region of Fig. 5. The radius of 1.27 ± O.O5-R 
from transit observations by Jha et al. (2000) is slightly outside, and the value 
of 1.146 ±O.O5Oi? from HST observations (Brown et al.) is slightly inside our 
3a confidence region (both times assuming a stellar mass of 1.1M Q ), with error 
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Table 2 

Determinations of planet and star sizes in the HD 209458 system 



Reference 


-Rpi(-Rjup) 


R*(R Q ) 


Rp\/R* 


Charbonneau et al., 2000 


1.27 ±0.02 


1.10 


0.1193 ±0.002 


Henry et al., 2000 


1.42 ±0.10 


1.15 


0.1276 ±0.009 


Mazeh et al., 2000 


1.40 ±0.17 


1.20 


0.1205 ±0.015 


Jha et al., 2000 


1.55 ±0.10 


1.27 


0.1261 ±0.008 


Brown et al., 2001 


1.347 ±0.06 


1.146 


0.1214 ±0.002 


this work 


1.435 ±0.05 


1.20 


0.1235 ±0.004 



bars from the HST observations well overlapping with our determination. 



3.3 Planet size, transit latitude and inclination 

In Fig. 4, the confidence regions along R v \/ and R'^ are shown. It should 
be noted, that the relative planet size R v \/R* is a function of R' # , whereas 
the absolute size R v \ is not. The region of fits with any significant meaning 
(er < 3cr m i n ) is limited to a range of stellar and planetary radii given by: 
R't = 1.20±0.06-Rq and R p \/R* = 0.1235±0.004. For a stellar radius of l.2R Q , 
the derived planet size would then be 1.435 ± 0.05i?j up . There is considerable 
scatter in the literature among derivations for the absolute size of the planet 
(Table 2), which is mainly caused by differences in the assumed stellar size. 
The agreement is much better if relative planet sizes are considered, whose 
average from all listed publications is R v \/ R* = 0.1231 ± 0.003, which is very 
close to the value derived in this work. 

The best fitting planet sizes in each of the four Stromgren colors are shown 
in Fig. 6. Their scatter around the mean is better than ± 3 %. This shows 
that the depths of the transit across the different colors agree with each other, 
and that the magnitude zero points in the light curves were correctly set. The 
remaining differences among the individual light curves will then primarily be 
caused by the dependence of limb-darkening on the wavelength. 

Independently of the stellar mass, the range of possible transit latitudes is 
restricted to b = 36.2 ± 6° (Fig. 4, within the 3 a confidence region). For 
= 1.1M , this corresponds to an inclination of the planet's orbit of % — 
85.96 ±0.55°, given by: 

cos? = — sin b (7) 
a p i 
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Fig. 6. Best fitting planet sizes in the four Stromgren colors, as a function of the 
stellar size (u: solid line; v. dotted; b: dashed; y: dot-dashed) 

where the planetary orbit's half axis a^M*, M pl , P) is given by Kepler's Third 
Law (a p i = 10.057i? o for M* = 1.1M ). 




4 The limb-darkening of HD 209458 



The excellent fit between the modeled light curves and data (Fig. 1) shows 
that the use of a linear limb-darkening expression has been adequate for this 
star. 

Values of the limb-darkening are shown for the four Stromgren colors in Fig. 7, 
again quoting w for the best fit at each value of R'^, similar to the dotted line 
in Fig. 3. Fig. 3 also shows that the derived values of w are not very dependent 
on small deviations of R p \ around its best value. Individual errors of the limb- 
darkening coefficients are estimated at ±0.05 for w v and Wb, ±0.07 for w y and 
±0.1 for w u . The derived values agree reasonably well with the usual sequence 
of limb-darkening coefficients along the four colors, with a distance between w u 
and w v that is smaller than their uncertainty. It is also apparent, that only for 
a small range of stellar radii, reasonable limb-darkening values can be assigned 
to HD 209458. For two narrow wavelength ranges, centered around 590nm and 
618nm, Brown et al. (2001) provide quadratic limb-darkening values, which 
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Fig. 7. Linear limb-darkening coefficients for HD 209458 in the Stromgren colors 
(u: solid line; v. dotted; b: dashed; y: dot-dashed) as derived from the fits. They 
are plotted against the normalized stellar radius of Eq.(6) and correspond to traces 
along best fitting values of R p \/R*, similar to the dashed line in Fig. 3. 

correspond to linear coefficients of w 59 o ~ 0.62 and w e i 8 ~ 0.58. These wave- 
length ranges are slightly red- wards of the Stromgren y filter (A = 547nm), 
and their limb-darkening agrees very well w y ~ 0.62 from our data, assuming 
a stellar radius of 1.146 R Q as did Brown et al.. 

The limb-darkening coefficients have been compared with values derived from 
the combination of evolutionary stellar models (Claret, 1995) and stellar at- 
mosphere calculations based on the ATLAS code, using tabulations of w from 
Claret (1998a). All the models have solar abundances and surface gravity val- 
ues were converted directly to stellar radii. Comparison with limb-darkening 
coefficients derived from PHOENIX code (Claret, 1998b) was not possible 
because the code does not cover the temperature range of HD 209458. The 
results from the evolutionary ATLAS models are shown in Fig. 8. Although 
a good general agreement between the limb-darkening sequence from transits 
and from ATLAS is found, the values from the transit have a wider spread 
across the four Stromgren colors. This can be seen through the intersections 
among the curves of w, which do not correspond to exactly the same stellar 
radii across the four colors. However, the intersections appear close to those 
stellar radii at which the best fits have been obtained. These intersection 
radii, averaged among the colors, are displayed in Fig. 5 (crosses) and have 
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Fig. 8. Limb darkening coefficients derived from the transit data (inclined curves) 
against values based on the ATLAS model atmosphere combined with evolutionary 
stellar models (Claret, 1998a, flat curves). Plots are shown for stellar masses of 1.05 
- 1.25 Mq, with the normalized radii (R#) of Fig. 7 converted to R* using Eq.(6). 
Stromgren filters correspond to line-styles as in the previous figures. 



a somewhat steeper dependency on M* than the radii directly derived from 
transits. For > 1.25M , no stellar model delivers radii small enough for 
those derived from transits, and for masses less than 1.05 M Q , radii from AT- 
LAS correspond to transit-radii outside of the 2a confidence region. The radii 
based on limb-darkening provide thereby a constraint for the stellar mass to 
be within 1.05 and 1.25 M , with a stellar mass of 1.1M leading to a radius 
of 1.17 ±0.03^0. 

The standard linear limb-darkening expression appears to be sufficient to ac- 
count for our observations. For one, HD 209458 is in the temperature range 
where the linear expression is expected to give good results. Second, most 
of the contributing non-linearity would come from angles close to the stellar 
rim, which are within the ingress and egress zone of the transit. Parameters 
describing non-linearities are therefore susceptible to introduce ambiguities 
with other transit parameters. Lastly, the number of free parameters in fits 
should be kept small in order to avoid the derivation of highly interdependent 
parameters of low significance. For stars outside of the temperature range of 
5000-6000 K, single-parametric limb-darkening expressions different to the one 
used here may be devised from modeled surface brightness distributions. An 



13 



Table 3 

Derived parameters for the HD 209458 system, and their interdependencies 



Parameter 


Value 


Assumption 


Reference 


R* 


1.20 ± O.O2i? 


M* = 1.1M 


Fig. 5, Eq.(5) 


Rpi/ R* 


0.1235 ±0.001 


Ri = 1.20 ± O.O2i? 


Figs. 4, 6 


Rpi/Rjup 


1.435 ± 0.05 


M* = l.lM ,i2* = 1.2O±O.O2i? 


Figs. 4, 6 


b 


36.2 ±2.2° 


Rp\/R* and R* from above 


Eq.(2), Fig. 4 


i 


85.96 ± 0.22° 


b and M* from above 


Eq.(7) 


w u 


0.82 ±0.10 


K = 1-1712© 


Fig. 7 


w v 


0.84 ±0.05 




a 


w b 


0.69 ±0.05 




a 


Wy 


0.59 ±0.07 




a 



example that this can be done with a very high precision, describing limb- 
darkening of the Sun, is given by Hestroffer and Magnan (1998). 

A possible planetary atmosphere may also cause differential color effects to be- 
come apparent in transit light curves. A color dependency of the atmospheric 
transparency, albedo, or scale height will cause the planet to display an ef- 
fective radius that changes with colors. We however do not expect to detect 
such color effects within the precision of our data. Assuming a precision of 
0.5 mmag, then SL ks 5A p i/A*, where A* = irRl is the cross section of the 
star, and 5A P \ is the differential cross section of the planet. With SA P \ depend- 
ing on a variation of the planetary radius like: 8A P \ = 2irR p \5R p i, we may solve 
for 5R pi as follows: 5R pi = 5A pl /2irR pl = R 2 JL/2R pl . For R pl = 1.43-R Jup = 
O.U8R Q and R* = 1.20R Q we obtain SR pl = O.OO241? w 1700A;m, meaning 
that color variations induced in the effective radius need to be larger than 
that value to be detectable in our data. Unless HD 209458 has an extremely 
extended atmosphere with a very unlikely transparency dependence along the 
Stromgren colors, our data will not contain any signatures from the planetary 
atmosphere. Data with significantly higher precision may however very well 
indicate the presence of planetary atmospheres. 



5 Conclusions 

A major goal of this work has been the derivation of limb-darkening coefficients 
for the star HD 209458. These, and other parameters of the system, have 
been obtained in a way which is quite independent of the underlying physical 
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assumptions. An overview of the derived parameters and assumptions is given 
in Table 3. References in the table are given to understand the interdependence 
of the derived parameters upon the assumed values. Any improvement in any 
parameter could be easily used, through the references given in the table, to 
adjust any other parameter for the HD 209458 system. For example, if we 
assume the stellar size of 1.146 R® from HST observations by Brown et al., 
and a mass of 1.1M , the corresponding planet size, transit latitude, and 
orbital inclination derived from our data would be R p \ = 1.323i?j up , b = 30.0° 
and i = 86.73°, which is in very good agreement with values cited by these 
authors. 

The derived limb-darkening coefficients are in good agreement with the AT- 
LAS atmospheric models, combined with evolutionary stellar models, over a 
stellar masses range within 1.05 to 1.25 M Q . There is a minor difference be- 
tween the stellar radii which best fit the transit light curve (1.20 ±0.02i? Q for 
M* = 1.1M ), and the radii which produce the limb-darkenings best describ- 
ing the ATLAS models (1.17±O.O3i? ). As the differences are within the error 
bars, we can not distinguish if they are caused by the errors inherent in the 
observational data, or by the calibration of the ATLAS models. More transit 
observations of HD 209458 with the same instrument are planned for the next 
year. It will then be possible to co-add light curves from individual transits, 
giving rise to an improvement in the precision of the parameters derived in 
this paper, especially with regard to the limb-darkening. 

With the discovery of further transiting systems, and with the availability of 
very high precision data from space missions, we expect that transit observa- 
tion will become an important tool in the verification of model atmospheres, 
and will be able to constrain stellar structure models over a wide range of 
stellar types. 
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